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The many-body states in an extended Fermionic Molecular Dynamics approach are flexible 
^ ■ enough to allow the description of nuclei with shell model nature as well as nuclei with cluster 
and halo structures. Diff'erent many-body configurations are obtained by minimizing the energy 
> \ under constraints on collective variables like radius, dipole, quadrupole and octupole deforma- 
^ \ tions. In the sense of the Generator Coordinate Method we perform variation after projection 
and multiconfiguration calculations. The same eff'ective interaction derived from realistic inter- 
^ ■ actions by means of the Unitary Correlation Operator Method is used for all nuclei. Aspects 
CO ■ of the shell model and cluster nature of the ground and excited states of ^^C are discussed. To 
understand energies and radii of neutron-rich He isotopes the soft-dipole mode is found to be 
important. 
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g : 1. FERMIONIC MOLECULAR DYNAMICS 
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In the Fermionic Molecular Dynamics (FMD) model [[H the A-body state is given as a Slater 
determinant | 2 ) of single-particle states | qi ) 

|g> = ^{|^i>^...^|^a>}. (1) 

The single-particle wave functions are described by Gaussian wave packets localized in phase- 
space 



<x|^> = J]c.exp{-^^^}|;^.>^ 



(2) 



In contrast to the AMD approach [ ^ in the FMD model the width parameter a is treated as 
a complex variational parameter that can be diff'erent for each single-particle state. The spin 
orientation given by the spinor \x) also treated as a variational parameter. For light nuclei 
the description can be improved by using a superposition of two Gaussian wave packets. 
The FMD many -particle state is determined by minimizing the intrinsic energy 

T (eK^Z:^ (3) 



with respect to the parameters of all single-particle states. This corresponds to a Hartree-Fock 
calculation using the FMD single-particle basis. To improve upon this mean-field result we 
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restore the symmetries of the Hamiltonian and perform a parity and angular momentum projec- 
tion. This defines the projection after variation (PAV) result. 

The eff'ects of the projection can be quite large and a variation after projection calculation 
(VAP) should be done [ 2|. But as a full VAP calculation is very expensive we perform VAP 
calculations in the spirit of the Generator Coordinate Method. The intrinsic energy is mini- 
mized under constraints on collective coordinates like radius, dipole, quadrupole or octupole 
moments obtaining a set of many-body configurations. The VAP minimum is then obtained 
by minimizing the projected energies of these configurations as a function of the constraints. 
These configurations can then also be used in a multiconfiguration calculation (MC) where the 
projected Hamiltonian is diagonalized within the set of these configurations. This provides us 
with a well defined procedure to study the properties of the ground and the excited states. 

2. EFFECTIVE INTERACTION 

We use an eff'ective interaction that is derived from the realistic Bonn or Argonne interac- 
tions by means of the Unitary Correlation Operator Method (UCOM) [ 3, 4J. The correlated 
interaction includes the short-range radial and tensor correlations induced by the repulsive core 
and the tensor force. The correlated interaction no longer connects to high momenta and can 
thus be used directly within product- state model spaces. The eff'ects of three-body correlations 
and genuine three-body forces are at this stage simulated by an additional two-body force that 
contains momentum-dependence and spin-orbit forces. This correction term is adjusted to re- 
produce the binding energies and radii of "^He, ^^O, "^^Ca, ^^O and "^^Ca. Here we find for ^^O and 
"^^Ca that after angular momentum projection the tetrahedral configurations (see Figure [T]) with 
Qf-clusters are favored energetically versus the spherical shell model configurations by about 
5 MeV Altogether we need a 15% correction to the ab-initio correlated two-body interaction. 




Figure 1. The eff'ective interaction is fitted to tetrahedral ^^O and "^^Ca and spherical ^^O and 
"^^Ca configurations. 



3. GROUND AND EXCITED STATES OF ^^C 

The question of the clustering in ^^C has a long history and has been studied in the AMD 
framework [|5l|6|. Our eff'ective interaction is fitted to binding energies and radii of "^He and ^^O 
but has not been adjusted to ^^C. We perform the FMD variation and obtain the (I s 1/2)^ (I 12)^^ 
shell model configuration (V), displayed in the left part of Figure |2l As it is spherical the an- 
gular momentum projection (PAV) does not aff'ect the result. The binding energy obtained with 
this configuration is 7.5 MeV too low compared to the experimental value (see Table [T]). The 
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obtained radius is too small. We can also restrict the intrinsic state to pure a-cluster config- 
urations. In these configurations the spin-orbit interaction does not contribute and already in 
the mean-field calculation a triangular configuration is preferred. In a VAP calculation where 
the positions and the widths of the a's are the generator coordinates we obtain 4.3 MeV less 
binding than with the shell model configuration and the radius is too large for the pure a-cluster 
configuration. If we perform a variation after projection calculation with the radius and the oc- 
tupole moment as generator coordinates we obtain the triangular configuration with polarized 
a-clusters denoted by VAP on the l.h.s. of Figure |2l It represents an intermediate between 
the VAP (2-cluster and the shell model (V/PAV) shape. The description of this polarization is 
strongly aff'ected by using two Gaussians per single-particle state. We obtain 4.2 MeV more 
binding compared to a VAP calculation with only one Gaussian per single-particle state. In 
case of the shell-model configuration the energy gain is only 1.5 MeV. The radius obtained in 
the VAP calculation is a little bit too small when compared to experiment. For the intrinsically 
deformed VAP configurations we can project out not only a 0"^ but also 2"^, 4"^ and 3" states (as 



Table 1 

Binding energies, charge radii and 5£'(2)-values of ^^C. 

^^[MeV] r,,,,,, [fm] B(E2) [eHm'] 



V/PAV 


84.7 


2.33 




VAP Qf-cluster 


80.4 


2.66 


56.3 


VAP 


91.9 


2.38 


24.7 


Multiconfig(4) 


93.4 


2.50 


40.0 


Exp 


92.2 


2.47 


39.7 ± 3.3 
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shown in the VAP column of the spectrum in Figure |3) and can calculate the B(E2) value for 
the 0"^ 2^ transition using bare charges. The B(E2) value obtained for the VAP configuration 
is too small compared to the experimental value. 

To improve the description we perform multiconfiguration calculations. For the additional 
configurations quadrupole and octupole moments are used as generator coordinates. In Fig- 
ure|2lthe configurations giving the biggest contributions for the first 3" and the second and third 
0"^ states are shown. In a multiconfiguration calculation containing the VAP and these three in- 
trinsic states the binding energy of the ground state is increased by 1 .5 MeV For the electromag- 
netic properties the eff'ect of admixing these less compact configurations is more pronounced. 
The radius is increased and the B(E2) value is now in agreement with the experimental value. 
The spectrum obtained with these four intrinsic states is shown in the Multiconfig(4) column of 
Figure 121 Using 25 configurations the energies for the ground state and the lowest excited states 
change only slightly. For the higher excited states the eff'ect is much bigger. Comparing to the 
experimental spectrum our results for the excited 0"^ and 2"^ states are already quite reasonable. 
Whereas the second 0"^ state seems to be mainly built of three a configurations, the third 0"^ 
appears to be mainly ^Be plus a. A detailed study of these 0"^ and also 2"^ states will be done in 
the future. We expect that even more extended configurations will have to be considered, taking 
into account the proposed Bose condensed nature of these states around the three a threshold [ 

13. 

4. HELIUM ISOTOPES 

To calculate the properties of the He isotopes we performed FMD calculations with two 
Gaussians per single-particle state using the electric dipole and the mass quadrupole moment as 
constraints for variation after projection (VAP) and multiconfiguration (Multiconf) calculations. 
In Table El the calculated binding energies and radii are shown together with the experimental 
values. In FigurejSlthe intrinsic configurations corresponding to the VAP minima are displayed. 
With the two Gaussians the FMD wave function is able to describe the neutron halo around 
the a-core. Using only a single Gaussian an imperfect halo together with a distorted a-cort is 
observed. Compared to the calculation using a single Gaussian the neutron-rich He isotopes 



Table 2 



Experimental and calculated binding energies and radii for the He isotopes. 







^He 


^He 


^He 


^He 


«He 


PAY 


Eb [MeV] 


28.3 


25.4 


26.3 


26.5 


29.8 


VAP 


Et [MeV] 


28.3 


26.8 


27.7 


27.7 


31.1 


Multiconf 


Et [MeV] 


28.4 


27.4 


29.1 


28.7 


31.7 


Exp 


Et [MeV] 


28.3 


27.4 


29.3 


28.8 


31.4 


Multiconf 


rmat [fm] 


1.45 


2.74 


2.42 


2.62 


2.53 


Exp 


rmat [fm] 






2.48 + 0.03 




2.52 + 0.03 


Multiconf 


rch [fm] 


1.69 


2.10 


2.02 


2.07 


2.03 


Exp 


rch [fm] 


1.68 
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Figure 3. L.h.s.: intrinsic shapes corresponding to VAP minima. Cuts through nucleon density 
at z = 0. R.h.s.: proton, neutron and nucleon densities as a function of radius. 



gain about 3-5 MeV in the PAV case. For the ^He nucleus we find a quadrupole deformed 
VAP minimum where the two halo neutrons are sitting on opposite sides of the core and a more 
bound (by about 1.1 MeV) dipole deformed VAP minimum where the neutrons are located 
on the same side of the core (see middle row of Figure [3l). For the other isotopes we only 
find minima with respect to the dipole deformation. The multiconfiguration calculations that 
include shapes of diff'erent electric dipole moments reproduce the experimental staggering of 
the binding energies. The large measured matter radii of ^He and ^He are also well described 
by the multiconfiguration calculations. We also observe a significant increase in the charge radii 
of the neutron rich He isotopes compared to "^He. This is explained by the zero-point motion of 
the more or less unchanged a-cort with respect to the center of mass of the nucleus. 

We conclude that the zero-point motion of paired neutrons in a soft-dipole mode seen in VAP 
and multiconfiguration mixing calculations is an essential ingredient for the understanding of 
neutron-rich He isotopes and probably also for other neutron halo nuclei. 
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